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We analyze the CCFR collaboration iron target data on the xF^ structure function 
k>( I making particular emphasis on the extraction of the higher twist contributions from 

data. Corrections for nuclear effects are applied in order to extract data on the 
structure function of the isoscalar nucleon. Our analysis confirms the observation 
made earlier, that the higher twist terms depend strongly on the level to which QCD 
perturbation theory analysis is applied. We discuss the impact of nuclear effects on 
the higher twist term as well as on the QCD scale parameter A-p^ extracted from 
the fit to data. 
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1 Introduction 

In the present paper we report the results of our analysis of the CCFR collaboration 
data |l|] on the structure function XF3. The particular emphasis of the analysis is to 
constrain the higher twist (HT) contributions to the structure function from data. The HT 
effects in the xF^ structure function are of particular interest because of certain theoretical 
predictions made in the framework of infrared renormalon technique 0, H, P, 0. An attempt 
to constrain the HT terms from the CCFR/NuTeV collaboration data was done in [|, Q, 
where the F3 structure function was written as the sum of two terms, 

xF3(x,Q2)=xF3^^(x,g^) + ^, (1) 

with F^^ the leading twist contribution and h/Q"^ the HT term. An important observation 
which follows from the analysis P, P, |lOl is that the magnitude of the HT term depends on 
the level to which the perturbation theory analysis of F^'^ is applied. If F^'^ is evaluated 
in a leading order (LO) renormalization group formalism a large h{x) appears from the 
fit. When a next-to-leading order (NLO) formalism is used for F^'^ a somewhat smaller 
but still substantial contribution from the HT term is needed. If F^'^ is evaluated to 
next-to-next-to-leading order (NNLO) very little room is left for the HT term. 

We note that QCD analysis of data implies that data are given for isolated proton and 
neutron. In practice, due to the reason of statistics, neutrino data are taken mainly on nu- 
clear targets rather than on isolated proton and neutron. For example, the CCFR/NuTeV 
collaboration uses the iron target, the IHEP-JINR Neutrino Detector uses the aluminum 
target [^, and the forthcoming data from CHORUS collaboration is obtained on the lead 
target [^). It is known from muon and electron DIS experiments, that nuclear effects are 
quite essential in a wide kinematical region of x and Q^ (the EMC effect at large x, nuclear 
shadowing at small x, for a review see e.g. [|I2|). Therefore, the separation of nuclear effects 
from data introduces certain corrections to QCD analysis of data. 

All these motivate us to make a new analysis of the CCFR neutrino data taking into 
account corrections due to nuclear effects. Our analysis involves two steps. In section |^ 
we discuss our approach to calculate nuclear structure functions and to correct data for 
nuclear effects, and then in section ^ we report the results of QCD analysis of corrected 
data. In sect. ^ we summarize. 



2 Nuclear structure functions 

In order to apply corrections for nuclear effects in our analysis we first calculate the "EMC 
ratio" for the iron target, R^^x^Q"^) = F^{x,Q'^)/AF^{x,Q'^), with ^3^ the structure 
function of a heavy nucleus of A nucleons and F^ the structure function of an isolated 
isoscalar nucleon.Q Then we extract the structure function of an isolated isoscalar nucleon 
from the CCFR data, Fi^{x,Q^) = F^^^^{x,Q^)/R3{x,Q^). 

Bulk of neutrino data with Q^ > 1 GeV^ is located in the region of a; > 0.1. For this 
kinematical regime it is usually assumed that nuclear DIS of leptons from nuclear targets 
can be viewed as incoherent scattering from bound nucleons. Major nuclear effects found in 
this region are due to nuclear binding |]I3| and Fermi motion |]I^ and off-shell modification 



of bound nucleon structure functions [|l^]. For the simplest nuclear system, the deuteron, 
the relation between the deuteron and the nucleon F3 structure function reads as follows 



[0, 



xF,^{x,Q') = 2/ ^ \^n{p)f (1 + :^) x'F,^{x',Q';p% (2) 

where ^E'_d(p) is the deuteron wave function which describes the probability to find the 
bound proton (or neutron) with momentum p, x' = Q'^/2p-q is the Bjorken variable of the 
bound nucleon with the four-momentum p which is given by the difference of the target 
four- momentum and the four- momentum of the spectator nucleon. Eq.(0) is written for 
the target rest frame and the axis z is chosen along the direction of momentum transfer, 
q = {qo,0±,—\q\). In this reference frame p = {M^ — y/p'^ + M'^.p) with Mr, and M 
the deuteron and the nucleon mass respectively and 7 = \q\/qQ = (1 + 4x^M^/(5^)^^^ is 
the 'velocity' of the virtual boson. Note that the bound proton and neutron are off-mass- 
shell and their structure functions depend on the nucleon off-shellness p^ as an additional 
variable. 

For the scattering off a heavy nucleus of A nucleons, there appears a rich spectrum of 

spectator nuclear states of A—1 nucleons, over which we have to sum. The nuclear structure 

function is then given by equation similar to (^ where we have to substitute the deuteron 

wave function by nuclear spectral function V{e, p) and introduce an additional integration 

^ The isoscalar nucleon structure function is defined as F^ ~ \{Fz + ^3*)- 



over the energy spectrum of spectator states TO 



xF,^{x,Q') = E /l^^^(^'P) (l + ^)x'F^{x\Q';p'), (3) 



de (Pp 

where the sum is over protons (r = p) and neutrons {t = n). The nucleon four-momentum 
p = (M + £:,p). The proton and neutron spectral functions, V^ and P", are normahzed to 
the number of bound protons {Z) and neutrons (A^) respectively. 

Heavy nuclei, such as iron ^^Fe26, generally have got unequal numbers of protons and 
neutrons with an excess of the latter over the former. The neutron excess is generally small, 
{N—Z)/A <^ 1. Therefore, it is a good approximation to assume that the neutron and the 
proton spectral functions calculated per one particle are equal, V^/Z = V"'/N. Then we 
find from (Q), 



where the averaging is done with respect to the isoscalar spectral function, pp+P". The 
last term in (^ gives a correction due to excess of neutrons in a nucleus. We notice that the 
sign of this correction is different for neutrino and anti-neutrino scattering. Indeed, we have 
F!^ = 2{d — u) with d and u the parton distributions of corresponding quarks in the target 
(we neglect for simplicity the contributions due to s- and c-quarks). Since the neutron 
has more (i-quarks than the proton has (in the valence quark region), the neutrino-neutron 
structure function is larger than the proton one, F^"^ > F^^. Therefore F^^ receives a 
positive correction due to excess of neutrons. Repeating this argument for anti-neutrino 
scattering, we find that the corresponding correction is equal in magnitude but opposite 
in sign (i.e. negative). Therefore the N—Z correction vanishes for the structure function 
averaged over neutrino and antineutrino. Similar discussion can also be applied to the F2 
structure function. One can also see that the N—Z correction is negative for the charged 
leptons scattering, i.e. for the F2 structure function. 

As it is obvious from (^, calculation of nuclear structure functions requires the knowl- 
edge of nuclear spectral function. In the next section we discuss the nuclear spectral 
function used in the present calculation in more detail. 

2.1 Nuclear spectral function 

Nuclear spectral function V determines the probability to find the nucleon with the mo- 
mentum p and (non relativistic) energy e in the ground state of the nucleus and can be 



written as follows 

Vie,p) = 2nJ2\{iA-lU-p\aM\A)fs(^e + Et' + j^^-E^y (5) 

Here the sum is over the quantum numbers of the whole set of the residual states of A—1 
nucleons which includes the bound states as well as the sates in continuum, a^ip) is the 
annihilation operator of the nucleon with momentum p and polarization a, and E:^~^ and 
Eq are respectively the energy of the residual nucleus and the ground state energy of 
the target nucleus. The residual system balances momentum of the removed nucleon and 
acquires the recoil energy p^ /2Ma^i though its effect is small for heavy nuclei. The nuclear 
momentum distribution is 

/ds 
zvr 

The integration of the spectral function over energy and momentum gives the number of 
bound nucleons A. 

The spectral function (^ determines the rate of nucleon removal reactions such as 
A{e, e'p)X that makes it possible to extract the spectral function from experimental data.0 
The picture of the spectrum of residual states as revealed from these experiments with heavy 
nuclei can be summarized as follows. For small energies and momenta^ the energy spectrum 
of residual states follows to that predicted by the mean field model of the nucleus, i.e. it 
consists of the set of sharp peaks whose positions can be identified with the energies needed 
to separate bound nucleons from the occupied single particle levels in the nuclear mean 
field. The deviations from the mean field picture become significant at high momentum p 
and high nucleon removal energy e. The widths of the resonances increase as e increases 
as well as the positions of the peaks move from that predicted by the energy independent 
mean field model. At high energy the spectral function is dominated by contributions from 
the states with one and more nucleons in the continuum. These contributions are due to 
NN-correlations in the nuclear ground state and can not be accounted for within the mean 
field model. 



^ Though one should notice, that the direct connection between the cross sections and the spectral 
function holds only in the impulse approximation, and is destroyed by other effects such as final state 

interactions and meson exchange currents. 

•^ The momenta should be compared with Fermi momentum pp which is for heavy nuclei pp ~ 
300 MeV/c. The corresponding Fermi energy is of order ep ~ 40 MeV. 



2.1.1 Phenomenological model of spectral function 

The calculation of the nuclear spectral function for complex nuclei requires to solve many 
body problem. The latter is known to be a difficult task and presently can be done only 
within certain approximations. In our discussion we follow |T^ and consider a phenomeno- 
logical model for the spectral function which incorporates both the single particle nature 
of the spectrum at low energy as well as high energy and high momentum components due 
to NN-correlations in the ground state.0 To this end we separate the full spectral function 
d^) into two parts, 

V{e,p)=VQ{e,p)+Vi{e,p), (7) 

which correspond to contributions from low excitation energy intermediate states ("Pq) 
and high excitation energy states (Pi). The low energy part can be approximated by the 
sum of the energy (5-functions which pick the positions of the occupied single particle levels 
weighted with the corresponding wave functions squared. In practice we use an approximate 
expression instead where the sum over occupied levels is substituted by its average value, 

Po(^,p) = 27rno(p)5 (^e + E« + ^^^ , (8) 

with E^^' = E^~^ ~Eq the nucleon separation energy averaged over residual configurations 
of A—1 nucleons with low excitation energies, i.e. mean field configurations, and no(p) the 
corresponding part of the nucleon momentum distribution. 

The high energy part Vi is determined by excited states in (^) with one or more nucleons 



in the continuum. It was observed within many body calculations |[T8|, |T^] for a wide range 
of nuclei that nuclear momentum distributions at high momenta {\p\ > pp with pp the 
Fermi momentum) run parallel to the deuteron distribution udIp), 

nM^C'nDip), (9) 

where the normalization constant C^ incorporates the many body aspects of the problem. 
It was found |]T^ that the constants C^ increase from 2 for '^He to 4.5 for ^^Fe. Going 
to a larger mass number does not bring in more high momentum component, C^ = 5 for 
nuclear matter .0 



* We note that our definition of the spectral function is different from the one used in ||l^ , where the 
recoil energy was not included into the energy (5- function in (pf). 

^ One should note however, that the relation (0) does not hold at low momentum where ni{p) contributes 
only a little to the full momentum distribution. 



This observation finds a simple interpretation if one assumes that high momentum 
component is generated by ground state configurations with a correlated NN-pair with a 
small distance between the nucleons. One can expect therefore that the relative motion in 
the NN-pair is determined by the properties of the NN-interaction in the vacuum rather 
than by long range nuclear interactions, and the distribution in the relative momentum 
will be similar to momentum distribution in the deuteron. 

In terms of the spectral function Vi this corresponds to the assumption about the 
dominance of the contribution from the states with one nucleon in the continuum and the 
remaining A— 2 nucleons being in a state with low momentum and low excitation energy, 

\A-l,-p)^a\p,)\{A-2r,P2)S{Pi+P2 + P)- (10) 

The coprresponding matrix element in (^) is then determined by the wave function of the 
NN-pair embeded into nuclear environment, 

((A-2)*,P2 \a{p,)a{P)\A) = 4^rcl{k)^CMiPcM)KPl+P2+P)- (11) 

We assume here factorization into the wave functions describing the relative motion in 
the NN pair, il),.e\{k), with relative momentum k = {p ~ Pi)/2 and the center-of-mass 
(CM) motion of the pair in the field of A— 2 nucleons, ipcuiPcM) with p^-^ = Pi + P- In 
general ipcM depends on the quantum numbers of the state of A— 2 nucleons, however the 
corresponding dependence of the V'rei is weak. 

We substitute ( pUf) into (H) and sum over the spectrum of A— 2 nucleons states and 
obtain an approximate expression for Vi, 

^i(^,P) = C^T^) J d^Pid^PcMnre\{k)ncM{PcM)^{Pi+P~PcM) (12) 

\ 2M 2Ma-2 J 

Here riy-ei and ricu are the relative and the CM momentum distributions respectively and 
ii;(2) _ ^-4-2 _ jj^A jg ^]^g energy needed to separate two nucleons from the ground state 
averaged over configurations of A— 2 nucleons with low excitation energy. Note that the 
minimum two nucleon separation energy E^'^^ = Eq~'^ — Eq is of order of 20 MeV for 
medium range nuclei like ^^Fe. 

The factorization of the matrix element (|TT|) into the relative and the CM motion wave 
functions is justified if relative momentum in the NN-pair is large relative to the CM 



momentum of the pair. This can be written as |p| ^ IPcmI- This condition allows us to 
approximate (^2|) by taking the relative momentum distribution out of the integral over 
the CM momentum at the point k = p. Then we have, 

P.,..P) ^ ,2.K„,,„(,(.^(^.^^B..))^^^. (13, 

where the averaging is done with respect to the CM motion of the pair. From the latter 
equation it is clear that the high momentum part of nuclear momentum distribution is 
given by the relative momentum distribution in the correlated NN pair embedded into 
nuclear environment, ni{p) = nrei(p)- 

The characteristic momentum for the CM motion of the NN-pair is similar to the one 
in the mean field model. In fact the averaged CM momentum squared of the pair can be 
estimated from the balance of the overall nucleus momentum [|^, {(J^PiY) = 0, where the 
sum is taken over all bound nucleons and the averaging is performed with respect to the 
intrinsic wave function of the nucleus. This gives (Pcm) = 2(^4— 2)(p^)/(A— 1), with (p^) 
the mean value of the squared single nucleon momentum. Since, by our assumption, the 
CM distribution does not include high- momentum component, we should also exclude the 
contribution of high-momentum part in estimating (p^). We follow |17] and parameterize 
the CM momentum distribution of the correlated NN pair in the field of other A— 2 nucleons 
by a Gaussian distribution, 

^cm(Pca/) = {(^/-^f^^expi-aplM), (14) 

with the parameter a determined from the averaged CM momentum of the pair, a = 

3/(2(p^m))- 

Using ([I^) we find that the integration over the CM momentum in (|l^) can be done 

analytically and finally the result reads, 

Vi{e,p) = ni(p) — v/^(exp(-ap^i„) - exp(-apL^)) , (15) 

where pmm and pmm are respectively the minimum and the maximum CM momenta allowed 
by the energy-momentum conservation in ([l^) for the given e and p, 



2 M-2, , \^ 2 M-2, , \^ 

With pr=(^(-2Mie + E^^^)-^)Y\ 



We notice that px has the interpretation of the maximal allowed CM momentum in the 
correlated NN-pair in the direction transverse to p for the fixed e and |p|. Note that the 
separation energy e is negative, as it follows from its definition in (^. The condition p^ = 
determines the threshold value of e for the fixed |p|. 



In numerical evaluations we use the parameterizations for nQ{p) and ^i(p) of [p!7| , 
which fit nicely the results of many body calculation of nuclear momentum distribution. 
It follows from this calculation that low momentum part incorporates about 80% of the 
total normalization of the spectral function while the other 20% are taken by the high 
momentum part. The mean kinetic energy obtained from integration of the full momentum 
distribution Uq + rii for the iron nucleus is (p^)/2M = 31MeV (the share of the high- 
momentum component rii is about 20MeV). The two parameters, E^^^ and E^'^\ determine 
characteristic range of nucleon separation energy. We set E^"^"* = Eq~'^ — Eq = 20 MeV, 
and therefore neglect possible contributions due to excited states of A— 2 nucleons in (p!2D.p| 
In order to fix the parameter E^^^ we employ the Koltun sum rule [|T9], which gives the 
relation between mean separation energy (e), mean kinetic energy (p^)/2M, and the ground 
state energy per nucleon Eq/A. For the mean kinetic energy of 31 MeV the sum rule gives 
(e) ~ —50 MeV. By integrating our model spectral function we find the value E^^^ = 27 MeV 
which satisfies the Koltun sum rule. 

2.2 The EMC ratios R2 and R^ 

In Fig.|I| we compare the iron/deuterium ratios for the charged lepton structure function 
E2 and the neutrino and antineutrino averaged structure function F3, calculated with 



the model nuclear spectral function discussed above. Also shown are the BCDMS |21 



and the SLAG [^ data on the iron/deuterium F2 structure function ratios. In numerical 



calculations we use the CTEQ4 parameterizations for the nucleon parton distributions pO |. 

We see that the behavior of the ratio R^ is very similar to that of the ratio R2 (for large x 

and large Q^ this does not come as a big surprise, since both F2 and xF^ are determined 

by valence quarks in this region). A small difference between the R2 and R^ curves is due 

to the neutron excess correction. As it was discussed in sect. |1|, the F2 structure function 

^ The effect of A— 2 excited states would lead to an overall increase of nucleon separation energy. We 
believe, however, that concrete estimates of this effect would require us to go beyond the model discussed 
in the present paper 



receives a negative N—Z correction, while similar correction cancels out in the neutrino and 
antineutrino averaged structure functions. It is well known, that the depletion of nuclear 



structure functions at x < 0.7 is due to nuclear binding effect [|13|, while the rise of the 
ratios at large x > 0.7 is due to nuclear momentum distribution effect (Fermi motion). 

We recall also that bound nucleons are off-mass-shell. Off-shell effects in the structure 
functions appear as the dependence on the target invariant mass p^. Target mass corrections 
can be of two different kinds. First of all, we have to take into account 'kinematical' target 
mass dependence due to finite ]? jQ^ ratio. To this end we use the Nachtmann scaling 
variable [^ ^ = 2x'/(l + (1 + 4x'^p^/(5^))^^^ instead of the Bjorken variable x' . Other 
('dynamical') sources of p^-dependence of structure functions are also possible. In this 
respect we refer to a model where p^-dependence of structure functions appears in the 



leading order |15, 16 1. We note also here, that we take into account off-shell effects in the 



bound nucleon structure function in a way that it does not affect the number of valence 



quarks in the nucleon |T6[. The off-shell effect acts coherently with nuclear binding effect 
and leads to an additional suppression of nuclear structure functions at intermediate range 
of X. 

The ratio R2 follows quite closely to data on the EMC effect in the iron nucleus (see 
Fig. |1[). This gives us the confidence in our method to calculate the EMC effect in the F3 
structure function. 

3 QCD analysis and fit 

Our QCD fit proceeds as follows. The nucleon structure function F^{x, Q"^) is written as 
a sum of the leading twist and the high twist terms, (|I|). We parametrize xF^^ at some 
scale Q^ = Ql in terms of a simple function, 

xF^'^i.x.Ql) =aix''^{l-xY^{l + aix). (17) 



Then we apply the renormalization group equation in order to calculate evolution of xFj 



LT 
3 

with Q^. We solve the renormalization group equation in the leading (LO), next-to-leading 
(NLO) and next-to-next-to-leadig (NNLO) logarithm approximations of QCD. In doing so 
we expand the leading twist structure function xF^'^ in terms of its Mellin moments within 
the framework of the Jacobi polinomial method and then apply the evolution equations to 



the moments (for more detail on the method used see |23|) 

10 



It should be noticed that in general the higher twist terms can be of two kinds: those 
which have the kinematical nature, e.g. the terms due to finite target mass, and those 
which arise due to higher twist operators and refiect the quark-gluon interaction effects 
in the target ("pure" higher twists). In order to ensure that the higher twist term in (|l]) 
describes effects due to quark-gluon interaction in the target we explicitly take into account 
the kinematical corrections due to finite target mass. To this end we substitute the Mellin 
moments by the Nachtmann moments [^ in the the Jacobi polinomial expansion of the 
leading twist structure function xF^^ . 

The CCFR data points are given in terms of discrete x-bins structure which range 
from X = 0.0075 to x = 0.75. We fit 116 data points with Q^ in the range between 
1.3 GeV^ and 200 GeV^. The fit parameters are the parameters a2, eta, and a^ of (|T^ at 
the scale Qq, the values of the function h{xi) at the center of each Xj-bin, as well as the 
QCD scale parameter A^p^. We fix the parameter ai by normalizing (p^Tj) to the Gross- 



Llewellyn-Smith sum rule, which was calculated in QCD to the second order in as ^5 
Sgls = 3(1 - «s/vr - 3.25(«s/vr)2). 

Our results are shown in Fig. ^ for the LO, NLO and NNLO approximations to the 
evolution equation. Also shown are the results with and without applying corrections for 



nuclear effects. We found that the fitting parameters are stable for Ql > 15 GeV^ and have 

?ll l_(2i\/^ Tr\v -i~rici T*CiCinl-fci T-\T*CiCiCi'n-fci*H ivi Hire r?l I ri cj t-wc 



chosen Qq = 20 GeV^ for the results presented in Fig. 0. The present fit includes more 



experimental points than that of |^, ^ . In particular, the inclusion of low-Q^ data points 
into the fit allows us to reduce the error bars in h{x) as compared to those presented in 
P, 1^. Though we should notice some increasing theoretical uncertainties associated with 
low Q^ data included into our analysis. We found that, in general, our present fit agrees 
with [^, ^], though introduces certain corrections especially at large x. 

A special care was taken to insure that our method to separate target mass correction is 
self-consistent. In particular we have done a special fit with h{xi) fixed at values presented 
in Fig. ^, but let the mass parameter in the Nachtmann moments to be free. We found 
that the minimum of x^ corresponds to the value of the mass parameter about 0.9 GeV 
with an error about 0.2 GeV. This value only weakly depends on the order of perturbative 
analysis and is close to the proton mass, that gives us confidence in the method used. 

As one can clearly see from Fig. |^ the magnitude of the higher twist effects depends on 
the level to which the perturbation theory analysis of F^^ is done. The more perturbative 
corrections are included into the evolution equation, the less room is left for the function 

11 



h{x). The shape of h{x) in NLO is in a quahtative agreement with the prediction of infrared 
renormalon approach |^. In particular, we found that the function h{x) is negative in the 
region 0.1 < x < 0.6. 

The separation of nuclear effects from data leads to a further suppression of the higher 
twist term h{x) at all levels of perturbation theory analysis of F^^. The effect of nuclear 
corrections on h{x) is most pronounced at large x, where we observe a systematic reduction 
of h{x) as compared with no-nuclear-effects analysis. Nuclear corrections result in the 
decrease of the values of the function h{x) at x > 0.6. As one can see from Fig. 0, the 
central points of h{x) at x = 0.65 and x = 0.75 bins become negative in contrast to the 
infrared renormalon prediction for large x [^]. An attempt to take into account nuclear 
effects in the QCD fit was previously done in |TT|, M. We comment in this respect that 



authors of |]Tl| used the deuteron model for nuclear effects, which is not a realistic one 



for the iron target. In particular, we found that h{x) is negative for large x, while it was 



positive in |jTT[. Authors of |^ attempted to introduce nuclear corrections to QCD fit in 
terms of the moments of structure functions. However, it was incorrectly assumed in 0, 
that the nuclear structure function F^ ^ as x ^ 1, that, to our mind, caused the x^ 
increase in their QCD fit. 

We found that the scale parameter A;^g^ decreases for about 40 MeV after nuclear effects 
are taken into account. This will lead to a shift of as{Mz) for about 2 ■ 10~^. Within the 
NNLO fit we get A^g^ = (394 ± 55) MeV for four quark flavors. 

4 Summary and conclusions 

In the present paper we report the results of our QCD analysis of CCFR data on F3 struc- 
ture function. The main emphasis was put on the extraction of higher twist contribution 
from data. We took special care to separate nucler effects from data, and compare the 
results of both analyses with and without corrections for nuclear effects. 

We found that nuclear effects cause about 10% decrease in the Ajj^ value. 

Our analysis confirms the observation made earlier, that the magnitude of higher twist 
terms decreases strongly when going from LO to NLO, and then to NNLO, approximations 
to the evolution equation. We observe an additional suppression of higher twist terms when 
corrections due to nuclear effects have been applied. 

In conclusion we note that small-x region in F3 structure function is of particular inter- 

12 



est, where a strong nuclear shadowing effect is anticipated PB|. We plan to address nuclear 



shadowing effect in application to QCD analysis of neutrino data. 
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Figure 1: The iron/deuterium ratios (EMC ratios) calculated for the structure functions 
F2 (solid curve) and F3 (dashed curve) within the model described in the text. The data 
points are from BCDMS |^ (open boxes) and SLAC |]22| (filled boxes) experiments. The 
curves were calculated at fixed Q^ = 16 GeV^ using CTEQ4 parameterizations for the 
nucleon parton distributions. 
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Figure 2: The function h{x), which describes the strength of the higher twist term in the 
X-F3 structure function as extracted from the fit to the CCFR neutrino data (see text). The 
labels on the figure indicate the level to which the perturbation theory analysis of xF^^ is 
done. 
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